Finite Density of States in a Mixed State of d x 2_ y 2 + id xy Superconductor 
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We have calculated the density of states of quasiparticles in a d x 2_ y 2 + id xy superconductor, and 
show that in the mixed state the quasiparticle spectrum remains gapless because of the Doppler 
shift by superflow. It was found that if the d xy order gap Ai oc yH as suggested by experiments, 
then thermal conductivity re oc \f~tL in accord with experimental data at lowest temperatures. This 
is an appended version of the paper published in Phys. Rev. B 59, 6024, (1999). We now also 
discuss the disorder effects and analyze the H log H crossover at small fields. We argue that H log 
H regime is present and disorder effect is dominant as the field-induced seconary gap is small at 
small fields. 
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Based on the experiments recently carried out by sev- 
eral groups , it has been found that the longitudinal 
thermal conductivity k xx of unconventional superconduc- 
tor, such as Bt2Sr2CaCu20g,j r s, displays quite strange 
behavior in the mixed state. With temperature from 5K 
to 20K, it decreases as the applied magnetic field along 
c-axis increases. At some critical field value Hk(T), ther- 
mal conductivity becomes insensitive to the increase of 
magnetic field and develops a plateau. In order to explain 
this phenomena, it has been proposed that the d x 2_ y 2 or- 
der parameter is suppressed at Hk due to the opening of 
a d x 2_ y 2 +id xy gap on the whole Fermi surface, thus sup- 
pressing quasiparticle contribution to thermal transport. 

On the other hand, more recent measurements at low- 
est temperatures say 0.1K, show thermal conductivity 
increases as field rises f^xx 

yH, as well as a substantial 
history dependence of the re in the mixed phase. These 
results clearly point towards the relevance of the vortex 
lattice and possible disorder and hysteresis in the vortex 
lattice for the thermal transport in the mixed state. 

From these measurements one can conclude that the 
spectrum of quasiparticles has a finite density of states 
N(E = 0, H) at lowest energies. While, we argue, based 
on the available data for re at lowest temperature, one still 
can not rule out the transition into fully gapped state, 
such as d x 2_ y 2 + id xy . 

Here we present the simple model calculation which 
indeed shows that for fully gapped state one can still get 
re oc \[H if Ai ex \fH as proposed by Laughlin [Q, and 
can reconcile the fully gapped spectrum of quasiparticle 
in d x 2_ y 2 +id xy state and the experimental data at lowest 
temperature by Aubin et al 0]. 

We repeated Volovik's calculation ||, in which the 
spectrum of quasiparticles is modified because of the 
Doppler shift from superflow in the mixed state. Volovik 
found that with the pure d x 2_ y 2 order parameter, the 
density of states of the delocalized quasiparticles has 
relation N(E — 0, H) oc yH, Here we consider 
d x 2_ y 2 + id xy order parameter, thus the quasiparticle 
spectrum is fully gapped, E(k) = ^/e 2 (k) + |A(k)| 2 , 



where A(k) = A (k) + iAi(k) , A (k) represents for the 
d x 2_ y 2 order parameter, and Ai(k) for d xy . Below we as- 
sume that there is some average value of the gap Ai (k) in 
the vortex lattice unit cell and ignore its position depen- 
dence. The spatial variations of Ai(k) albeit substantial, 
will not change qualitatively results presented below. 

The DOS for delocalized quasiparticles at Fermi sur- 
face can be calculated from following equation, taking 
T = 0, 
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where v s = (h/2m e )((f>/r). The main contribution of the 
DOS comes from the vicinity of the d x 2_ y 2 gap nodes, 
the gap function has the general form around the node 
A (k) <f(fc z ) • [k - k n (k z )], where 7(fc z )||e z x k„.Jt is 
convenient to introduce new momentum variables k x — 
(k — k„)-fc n , k y = (k— k„)-7„, where k n is the unit vector 
in the direction of the gap nodes in the a — b plane, j n 
the unit vector along the direction of 7. Then we can 
write Eq.(l) as 
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5[yfe 2 + k2~f 2 (k z ) + Ai 2 + m e v F (k„) • v s 
and the integral can be simplified as 
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Substituting \J x -V Ai 2 by the variable y, we can further 
write above equation as 
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In order to find the analytical form we first express the 
delta function as an integral form, S(y + kpv s cos9) = 
dXexp (iXy — i\kFV s cos9), and integrate over 9, 
then over A. We find from Eq.(4), 
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where ro = |£A _1 , A = Ai/Aq, Rh is the intervortex 
distance, and £ is the coherent length for pairing. We al- 
ready know that Rh ~ £\/ H c2 /H, and assume Ai much 
smaller compared with Ao, hence Min(ro, Rh) = Rh- 
After making Taylor expansion in terms of Rh/tq, we 
can get the following relation averaged over the vortices 



A 2 

y/W-^L. (6) 
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Here the factor K is on the order of unity. We found that 
the spectrum remains gapless. The origin of the gap- 
less behavior is the Doppler shift of quasiparticle states 
E — kv s , which is position dependent. There are regions 
where this shift is larger than minimal gap in the spec- 
trum Ai, thus leading to the finite DOS. This point was 
also emphasized by Hirschfeld and Wolfle, who showed 
that for finite id' gap the density of states will be finite 
H . We note that regardless of the power with which 
Ai ~ H a opens up in the field, as long as a > 1/2, the 
leading term in DOS will always be y/H at small fields. 

If one uses the form of Ai = twy/2eH/{hc) used by 
Laughlin [Q (v is the average quasiparticle velocity in the 
d-wave node) to explain Krishana's data, we find that for 
the DOS of the fully gapped state 
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where Ep is the Fermi energy. Typically A /Ep ~ 10 1 
and the effect of induced gap is to lower the y/H prefactor 
in N(0) by about 20%. The DOS with N{E = 0, H) ~ 
y/~H was previously thought to be characteristic of "pure" 
d x 2_ y 2 state ||. 

Experimental result by Aubin et al. showed that 
the thermal conductivity at lowest temperature with 
k(H) oc \f~H_ is consistent with quasiparticle transport 
and "Volovik effect" in the density of states. They also 
pointed out that even at lowest temperature data are 
consistent with the gapped d x 2_ y i +id xy phase provided 
Ax is small enough. 
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FIG. 1. Two regions of the phase diagram are shown. The 
line is the the boundary betwe en gapped and gapless states, 
given by ksTc. = 0.52hv\J ^£lL . In both phases the density of 
states N(0, H) oc y/H due to A xy oc y/~H. 

Fig.l demonstrates that the superconductor has two 
phases [0 as temperature and applied magnetic field 
change, one with "pure" d x i_ y i order parameter, the 
other with d x 2_ y 2 + id xy order parameter. According to 
our calculation, in both phases, the DOS of quasiparticle 
at E = is proportional to square root of applied field. 

We also note the implication of our calculation for the 
specific heat. Recent measurements of the specific heat 
C(T, H) in the mixed state of YBCO superconductors by 
Moler et.al. § indicate that the low temperature density 
of states indeed scales as N(E = 0,H) ~ y/H. This 
observation however, as we argue, is not inconsistent with 
the d + id' gap in the mixed state. 

In the paper by Volovik Q, it was mentioned that in 
the mixed state of the d-wave superconductor, the con- 
ventional s-wave pairing will be generated in the core of 
vortex. This might be another possibility to obtain fully 
gapped quasiparticle spectrum. However, it will not af- 
fect our result which comes from the region outside the 
core. 

The main result of this work is to show that even in the 
presence of fully gapped spectrum, the superfiow mod- 
ifies quasiparticle spectrum and keeps DOS gapless in 
the mixed state. Approximating thermal conductivity as 
proportional to the DOS at E = at lowest temperature, 
and assuming that the magnitude of Ai proportional to 
y/H in accordance with the work by Laughlin, one can 
still be consistent with experimental data at lowest tem- 
perature. It shows that we cannot rule out the possibility 
that a second superconducting phase appears in the mag- 
netic field at low temperatures. 

Note Added (March, 1999): In the paper by C. Kubert 
and P.J. Hirschfeld 0, the density of states was calcu- 
lated in a d x 2_ y 2 superconductor in the dirty limit, and 
they found a H\nH behavior at the lowest field instead 
of y/H . Here we also calculated the density of states with 
the induced gap Ai in the presence of impurity scatter- 
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ing, and found that the appearance of the induced gap Ai 
does not change the HlogH behavior in the lowest field. 
The details of calculation are described in the appendix. 
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1. Appendix A 

We assume the impurity gives an imaginary term 70 in 
the self energy of one particle Green function. Thus the 
density of states can be calculated by following equation: 

N(E = 0,H) = -jL- f d 3 k f d 2 vIm{G) 
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FIG. 2. Normalized density of states vs (H/H c2 ) 1/2 . The 
three line groups for each data sets correspond to different 
strength of Ai = xA (H/ H c2 ) 1/2 with x = 0, 0.5, 0.8 from top 
line to bottom. The presence of the field induced Ai oc H 1 ' 2 
does not change the Nh (0) behavior at lowest fields as the ef- 
fect of disorder and Ao dominates this region. Hence for dirty 

d x 2_ y 2 + id xy we expect the same H 1 / 2 — » HlnH crossover, 

7pas discussed before M. The data are from references [6] and 
[8]. 
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where Eh = \/ H/ H^Aq, Y is an upper cutoff of 
e 2 . + I A| 2 . It should be noticed that 70 is the imaginary 
self energy term from impurity, while 7 is the gradient of 
Ao with respect to k y 

In above equation, we can get the HlogH behavior 
when Ai equals to zero. We calculated the integral nu- 
merically, including both the effect of induced d xy gap 
and impurity, and show them in Fig. 2. In our calcula- 
tion we used the result, 7 = 0.61\/rAo, of Kubert and 
Hirschfeld @, with r = ni/nNp- It is found that the 
presence of an induced d xy gap still cannot change the 
behavior of the density of states as a function of applied 
field, however it decreases the density of states a little 
bit. 
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